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INTRODUCTION 

Preliminary results [1] of Japanese snowmelt runoff modelling were proposed at First Japan/US 
Snow and Evapotranspiration Workshop held in Tokyo on March, 1982. Aiming at the daily 
snowmelt runoff forecasting, the preliminary model was designed through multiple regression 
analyses based on the hydrological conception. However, it involved several operational defects as 
follows: 

i) containing three input variables of mean air temperature, precipitation and snow cover area at the 

day of which the inflow ought to be forecast. 

ii) insufficient concordance between the hydrological conception and the parameters which resulted 

from the multiple regression analyses. 

In the present study the model was, then, so modified that all the input variables are of the 
antecedent days. The inflow of the previous day is also taken into account. The new model was 
developed in use of the 1980 data in the Okutadami River Basin, chosen again as the test basin, and 
was adapted to the other years for its operative efficiency. 

The optimum contribution of the Satellite imagery to the snowmelt runoff forecasting is the 
extraction of the snow coverage. Japan is, unfortunately, covered with much cloud in snowmelt 
season. Only in a fortunate season, one can obtain only one or two useful imagery. Unlike US it is, 
therefore, impossible in Japan to draw the depletion curve of snow coverage based on the imageries. 
However, a few imageries obtained in the past were effectively utilized for the verification and 
modification of the depletion curve induced from the snow water equivalent distribution at maximum 
stage and the accumulated degree-days at one representative point selected in the basin. This point is 
equipped with data collection facilities and is expected to fully function as DCP. 

Together with the depletion curve, the relationship between the basin-wide daily snowmelt 
amount and the air temperature at the point above were exhibited in form of nomograph for the con- 
venience of the model user. 

The Japanese model proposed in the present study is, thus, the practical model of the snow melt 
runoff forecasting. 

TEST BASIN 

The Okutadami River basin (area 414 km 2 ) chosen as the test basin for this analysis is situated 
in Central Japan shown in Figure 1. The highest and lowest altitudes of the basin are 2346 m and 
782 m respectively. The altitude difference is, therefore, 1564 m. This basin is one of the heavies't 
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snowfall area in Japan; the maximum snow depth recorded in recent 20 years was 598cm at 
Okutadami Meteorological Observing. Station located at the dam site in the basin. The station runs 
the measurement of snow water equivalent and inflow into the reservoir as well as the general 
meteorological observation. Located slightly outside the basin is the Hinoemata meteorological 
observing station of air temperature and precipitation. Two stations above and the representative 
point (mentioned later), Maruyama, provide all the necessary input data of the model. Those loca- 
tions are shown in Figure 2. 

The snow depth measurements by helicopter at ten to twenty points scattered in the basin are 
carried out at the maximum stage of snow depth at about the end of every winter season by the 
Electric Power Development Company. The measurement points by helicopter are shown in 
Figure 3. 

CONCEPT OF SNOWMELT RUNOFF MODEL 

Two forms of multiple regression models were considered. First model is as follows; 


R„ 

= c + a,X n _j + a 2 X n _ 2 + a 3 X n _ 3 

(1) 


+ b,Y n _, + b 2 Y„_ 2 + b 3 Y n _ 3 


X n _i 

= (M n -i + Pn-i^'n-i 

(2) 

Y n -i 

= P»-iS"„-i 

(3) 


where 


Rn 

: daily direct inflow (m 3 /sec) into reservoir on the day, n 

C 

: constant 

,b 2 ,b 3 

: partial regression coefficient 

M n _i 

: mean snowmelt amount on the day n— i (m 3 *sec -1 *km -2 ) 

P„-i 

: mean rainfall amount on the day n-i (m 3 *sec -1 *km -2 ) 

S'„-i 

: area where daily mean air temperature is above 0°C within snow cover area 
on the day n— i (km 2 ) 

S"„-i 

: area where daily mean air temperature is above 0°C within non-snow cover 
area on the day n— i (km 2 ) 


The second model, in which the direct inflow of the previous day, R,,.,, is added to the 
variables of the first model, is as follows; 

R„ = c + dR,., + ajX,., + + a 3 X n _ 3 (4) 

+ + b 2 Y n _ 2 +~ b 3 Y n _ 3 
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where 


d : partial regression coefficient 

DEPENDENT VARIABLE 
Direct Inflow 

The dependent variable of the model is the daily direct inflow, R,,, into the reservoir. The base 
flow which should be subtracted from the daily inflow to obtain the direct inflow, R,,, is expressed 
by 


r, = r 0 (1.0+0.00124t) 


( 5 ) 


where 

r, : base flow after t days from the day when first discharge occurs due to snowmelting in the 
season (m 3 /sec) 

r 0 : inflow before the snowmelt season starts (m 3 /sec) 

The equation (5) was empirically determined by the inflow data recorded at the beginning and end of 
the several snowmelt seasons 

INDEPENDENT VARIABLES 

Snow Cover Area and Its Related Areas, S n _„ S' n _| and S" n _| 

The snow coverage information from Landsat MSS is scarcely obtainable in Japan. Then, the 
snow line during snowmelt season was considered to appear when snow water equivalent in a mesh 
(shown in Figure 2) became equal to the snowmelt amount which can be calculated by the degree- 
day method (mentioned later). In order to extract the snow line, the computation of the snow water 
equivalent in each mesh at the maximum stage is required. Then, the computation was made by 
following three methods; the topb'graphic method [3] in use of the extensive snow survey data in 
1980, the degree-day method [4] in effective use of two Landsat MSS data fortunately obtained in 
1979 and the comparison method [5] in use of the regional snow depth data obtained for several 
years with helicopter. Next, the date when all the snow deposits disappear was calculated for each 
mesh, applying the degree-day method. The snow cover area and its related areas, s n _j, s and 
s" n _j were then counted for every day during each snowmelt season. 

Total Snowmelt Amount, M n .|*S' n _| 

The degree-day method was adopted for the snowmelt calculation, expressed by 

M = KET (6) 
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where 


M : snowmelt amount (g/cm 2 ) 

ET : accumulated daily mean air temperature above 0°C (°C*day) 
k : degree-day factor (0.45; average value in preceding studies) 

The daily mean air temperature at the Okutadami Meteorological Observing Station was adopted 
as the standard for the estimation of the air temperature at all meshes assuming the lapse-rate as -0.6 
°C/100 m. 

Applying the equation (6), the calculation of the daily snowmelt amount was made for each 
mesh within the area, S ' n _., during each snowmelt season. By summing up them, the total snowmelt 
amount in whole basin, M n _ i »S' n _ i was then obtained. 

Total Rainfall Amount, P n _,*S' n _| and P n _j*S" n _i 

The Okutadami Station is the only one rainfall observation station within the basin running 
throughout snowmelt season. However, located at the lowest elevation and the northern edge of the 
basin, it does not give the mean rainfall of the basin. Then, the Hinoemata Station was adopted for 
this purpose, although it is located at slightly outside the basin. After careful meteorological inspec- 
tion, it was assumed that the precipitation at Hinoemata is equal to the rainfall amount in each mesh 
where the daily mean air temperature is above 0°C both in snow cover and non-snow cover area. 

The total rainfall amount, P^-S'^ and P n _j*S" n _j were then obtained. 

CONSTRUCTION OF SNOWMELT RUNOFF MODEL 

Now, all the input data were prepared for the determination of the constant and partial regres- 
sion coefficients of the equations (4) and (5). The multiple regression analyses were, then, carried 
out in use of data in 1980. The resultant constructed two models are presented in Table 1. 

The comparison between the observed and calculated inflow, exhibited in Figures 4 and 5, 
shows that the models would stand the practical use the snowmelt runoff forecasting. 

NOMOGRAPHS FOR PRACTICAL USE OF MODELS 

In course of the model construction, the input variable values were obtained through the com- 
plex calculation procedure relating to the mesh. Necessary for the convenience of the practical use of 
the models is an easy acquisition of the input variable data. Following nomographs were, then, 
devised. 

Nomograph (I) for Basin-Wide Snow Amount at Maximum Stage 

The Maruyama shown in Figure 2 was selected among several tens proposed sites in the basin 
as the best representative site, whose snow water equivalent has a linear correlation with the basin- 
wide snow amount throughout both snowfall and snowmelt seasons. Based on the six years data is 
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Figure 6 (nomograph (I)), which exhibits a relationship between the basin- wide snow amount and 
Maruyama’s snow water equivalent at the maximum stage. 

The Maruyama is situated at roughly mean altitude of the whole basin and is at present equipped 
with DCP. 

Nomograph (II) for Snow Cover Area 

The depletion curve of snow cover area with increasing degree-day at the Maruyama was 
obtained for four snowmelt seasons, and is exhibited in Figure 7. The curves seem to depend on the 
snow amount at the maximum stage. Then, the relation was investigated between the snow amount 
above and the Maruyama’s degree-day accumulated up to the first day when the bare ground appears 
somewhere in the basin, namely, S n _;/A<l. A is the basin area. A linear correlation was found as 
shown in Figure 8. Based on this correlation, the depletion curve was standardized. The resultant 
curve shown in Figure 9 is the nomograph (II). 

Nomograph (III) for Basin-Wide Snowmelt Amount 

A clear relation was found between the basin-wide daily snowmelt amount and the daily mean 
air temperature at the Maruyama through the medium of the ratio (S n _j/A). This relation shown in 
Figure 10 is the nomograph (III) for M n _;— S' n _;. 

Nomograph (IV, V) for Basin-Wide Rainfall 

In the models, two kinds of basin-wide rainfall amount are taken into account; one in the snow 
cover area (S nH ) and another in the non-snow cover area (A-S n _j). 

Regarding the snow cover area, the extraction of the area (S'^n-i) where the daily mean air 
temperature is above 0°C is necessary. For simplicity, the relation was obtained between the total 
rainfall amount (P n _j»S' n _j) and the Hinoemata precipitation (P n _ ; ) through the medium of the pro- 
duct of the Okutadami daily mean air temperature and the ratio (S n _j/A). The relation shown in 
Figure 11 is the nomograph (IV) for P n _j*S' n _j. 

As for the non-snow cover area, this area is usually above 0°C during a snowmelt season, so 
that S" n _in— i -5- A — S n _j. Exhibited in Figure 12 is the nomograph (V) for P n _ i »S" n _ i . 

VERIFICATION OF SNOWMELT RUNOFF MODELS 

In order to verify the models (Table 1 constructed for the 1980 data, the models were applied to 
the other two years, 1979 and 1982. All the input values were obtained in use of the nomographs 
(I)-(V). The results are shown in Figures 13-16. A good agreement can be seen except May 26-28, 
1979. A large difference in this period may be caused by heavy rainfall unusually in a snowmelt 
season which occurred on May 26 and 27. In Japan the rainy season starts at the beginning of June, 
and often heavy rainfall occurs. The fact above shows that the models should not be applied in case 
the heavy rainfall happened near the end of snowmelt season when the snow cover area shrinks up 
to the top of the high mountains. 
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CONCLUSION 


In the present study scarce Landsat MSS data obliged the limited but efficient utilization for the 
verification of the depletion curve of the snow cover area. In order to construct the models, the 
complex topographic mesh analyses were carried out. However, the data of final necessity become 
only the mean air temperature and precipitation for the daily use of the models. 

The runoff forecasting procedure is summarized in Figure 17. 
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Figure 1. Location of The Test Basin 
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Figure 3. Locations of Points Where Snow Depth and Its Water Equivalent 
were Measured by Electric Power Development Co. 
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Table 1 

Snowmelt Runoff Model 


w 

00 



MODEL-II 

Rn = 7.33 * 0.549 Rn-,+ 0.392 Xn-i - 0.1 1 0 Xn- 2 

- 0,0 5 7 Yn-i 

m.c.c. : 0.943 s.d. : 14.7 
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Figure 6. NOMOGRAPH(I) for Basin-Wide Snow Amount at Maximum 
Stage in Use of Snow Water Equivalent at Maruyama 
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Figure 7. Relationship Between Depletion of Snow Cover Area and Accumulated Degree Days at 
Maruyama (Parameter: Basin-Wide Snow Amount at Maximum Stage) 
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Figure 8. Relationship Between Basin-Wide Snow Amount and Accumulated 
Degree-Days at Maruyama Before Snow Cover Area Depletion Starts 
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SNOW COVER AREA ( S~/a ) 



Figure 9. Nomograph(ll) for Snow Cover Area. (Parameter: Basin-Wide Snow Amount at Maximum 

Stage.) 
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Figure 10. Nomograph(lll) for Basin-Wide Snow-Melt in Use of Daily Mean Air Temperature at 
Maruyama. (Parameter: Relative Snow Cover Area, S n _/A) 
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m.c.c. : 0.835 
s.d. : 17.0 
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Figure 13. Comparison Between Observed (Solid Line) and Calculated 
(Broken Line) Inflow in Use of MODEL-1 for 1979 Snowmelt Season 
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Figure 15. Comparison Between Observed (Solid Line) and Calculated 
(Broke Line) Inflow in Use of MODEL-II for 1979 Snowmelt Season 
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Figure 16. Comparison Between Observed (Solid Line) and Calculated 
(Broken Line) Inflow in Use of MODEL-II for 1982 Snowmelt Season 




Figure 17. Procedure of the Snowmelt Runoff Forecasting 













